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METHOD FOR FAULT DIAGNOSIS OF A TURBINE ENGINE 

BACKGROUND OF THE INVENTION 

[001] The present invention generally relates to performing fault diagnosis 
of turbine engines and, more specifically, to an automated method for providing 
such fault diagnosis. 

[002] The gas turbine propulsion engine is one of the most critical 
5 subsystems within an aircraft. A great deal of research has been carried out in 
the field of gas turbine propulsion engine fault diagnosis and prognosis. The 
process of diagnosis requires adequate operating data input which is not 
always available or provided. The mid-sized jet propulsion engine, in particular, 
may typically be operated with less-expensive data acquisition systems, 

10 because of cost constraints incurred in design and operation, and the quality of 
operating data output suffers as a result. For example, engine output data may 
be recorded only during certain pre-determined operating intervals of time. 
Accordingly, only a limited number of parameters related to the core engine 
performance are routinely recorded and available for analysis. Moreover, 

15 conventional diagnostic methods typically fail to take into account operating 
differences resulting from engine-to-engine manufacturing variations. 
[003] Many of the heuristics derived from an expert's knowledge and 
understanding of engine performance are not systematic in nature. 
Accordingly, the current practice in the field has mostly relied on manual 

20 trending of the observations as performed by the experts. Although the present 
state of the art discloses unconventional methods for enhancing manual 
observation and control, such as the method incorporating a self-organizing 
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map for controlling process events of a technical plant taught by U.S. Patent 
No. 6,314,413 issued to Otte, such procedures are not generally applicable to 
the diagnosis and prognosis methods required in the maintenance and 
operation of the gas turbine engine. In addition, conventional methods of 
5 engine diagnosis may not provide a confidence level or belief factor for a 
determination that the engine condition is faulty or normal. 
[004] As can be seen, there is a need for an improved method of 
performing fault diagnosis and prognosis of turbine engines. 

1 0 SUMMARY OF THE INVENTION 

[005] In accordance with the present invention, a method of analyzing a 
turbine engine comprises acquiring an engine operating parameter; calculating 
an engine residual value from the engine operating parameter; normalizing the 

15 engine residual value to yield a normalized engine residual; mapping the 
normalized engine residual as an input vector into an engine condition space 
having a plurality of clusters, each cluster representing either a normal engine 
condition or a faulty engine condition; and identifying a closest cluster within the 
engine condition space to determine whether the engine under analysis is 

20 normal or faulty. 

[006] In another aspect of the present invention, there is provided a 
computer readable medium having computer-executable instructions for 
performing a method, wherein the method comprises calculating an engine 
residual value from data generated from an engine model and from engine 

25 operating data collected in the field; normalizing the engine residual value to 
yield a normalized engine residual; mapping, via a self-organizing map, the 
normalized engine residual as an input vector into an engine condition space 
having a plurality of clusters, each cluster representing either a normal engine 
condition or a faulty engine condition; and identifying a closest cluster within the 
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engine condition space to determine whether the engine condition is normal or 
faulty. 

[007] In yet another aspect of the present invention, a method of 
analyzing a turbine engine to determine a normal or faulty condition comprises 
5 acquiring a plurality of engine operating parameters; calculating a 
corresponding plurality of engine residual values by comparing each of the 
engine operating parameters with standard engine characteristics obtained 
from an engine model; computing the mean and the standard deviation of each 
engine residual value obtained from a normally operating baseline engine; 

10 normalizing each engine residual value by normalizing the mean to zero and by 
normalizing the standard deviation to unity to yield a plurality of normalized 
engine residuals, the step of normalizing dynamically using normalization 
factors obtained from a parameter distribution of a normally-operating baseline 
engine; mapping the normalized engine residuals as input vectors into an 

15 engine condition space having a plurality of clusters, each cluster representing 
either a normal engine condition or a faulty engine condition; and identifying a 
closest cluster within the engine condition space, the closest cluster being 
closer to the input vectors than any other of the plurality of clusters, the closest 
cluster representation determining whether the engine under analysis is normal 

20 or faulty. 

[008] In still another aspect of the present invention, a method of 
analyzing a turbine engine to determine a normal or faulty condition comprises 
training a self-organizing map by inputting data from a plurality of reference 
turbine engines; acquiring a core speed reading, an exhaust gas temperature 
25 reading, and a fuel flow reading for the engine under analysis; calculating a 
core speed residual value, an exhaust gas temperature residual value, and a 
fuel flow residual value by comparing the core speed reading, the exhaust gas 
temperature reading, and the fuel flow reading with corresponding standard 
engine characteristics obtained from an engine model; computing the mean 
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and the standard deviation of each of the core speed residual value, the 
exhaust gas temperature residual value, and the fuel flow residual value which 
are obtained from a normally operating baseline engine; normalizing each of 
the core speed residual value, the exhaust gas temperature residual value, and 
5 the fuel flow residual value by normalizing the respective means to zero and by 
normalizing the respective standard deviations to unity to yield a normalized 
core speed residual, a normalized exhaust gas temperature residual, and a 
normalized fuel flow residual; mapping, via the self-organizing map, the 
normalized core speed residual, the normalized exhaust gas temperature 

10 residual, and the normalized fuel flow residual as respective input vectors into 
an engine condition space having a plurality of clusters, each cluster 
representing either a normal engine condition or a faulty engine condition; and 
identifying a closest cluster within the engine condition space, the closest 
cluster being closer to the input vectors than any other of the plurality of 

15 clusters, the closest cluster representation determining whether the engine 
under analysis is normal or faulty. 

[009] These and other features, aspects and advantages of the present 
invention will become better understood with reference to the following 
drawings, description and claims. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 

[010] Figures 1A and 1B jointly show a flow diagram schematically 
representing a method for mapping acquired engine operating parameters or 
25 data values into an engine condition space from which an engine diagnosis and 
a belief factor can be obtained, according to the instant invention; and 
[011] Figure 2 shows a computing module or device which can be used to 
implement the method of Figures 1 A and 1 B. 
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DETAILED DESCRIPTION OF THE INVENTION 

[012] The following detailed description is of the best currently 
contemplated modes of carrying out the invention. The description is not to be 
5 taken in a limiting sense, but is made merely for the purpose of illustrating the 
general principles of the invention, since the scope of the invention is best 
defined by the appended claims. 

[013] Using the disclosed method for fault diagnosis of a turbine engine, 
aircraft maintenance personnel, for example, can establish whether a turbine 

1 0 engine under analysis is operating normally or may have a faulty condition. The 
maintenance personnel can input parameter values, obtained in the process of 
routine maintenance or obtained from flight data recorded in the field, into a 
diagnostic computer, for example, and derive a diagnosis or prognosis output 
along with a confidence level indication, or belief factor, for the analysis. The 

15 failure modes addressed may include failures related to the turbine engine 
combustor liner, bleed system, and exhaust gas temperature sensors. 
[014] The diagnostic process uses a mapping technique to automatically 
correlate the inputted parameter values to a probable engine operating 
condition selected from a plurality of normal and faulty engine operating 

20 conditions previously "learned" by the mapping technique, to derive the engine 
condition diagnosis and provide the confidence level to the maintenance 
personnel. 

[015] The disclosed method for fault diagnosis further allows aircraft 
operators to perform preventive maintenance, by anticipating potential faulty 
25 conditions, and thus achieve a reduction in engine operating costs. 
Conventional diagnosis methods, in comparison, may need to rely on the 
intervention of a highly-trained operator having expert knowledge and 
understanding of the engine operating characteristics because of the limited 
operating data normally obtained during aircraft maintenance. In addition, 
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conventional diagnosis methods typically do not provide a prognosis and 
confidence indication for the engine. Without such reliable prognosis, there 
remains the potential for aircraft flight delays and cancellations, unscheduled 
engine removals, and takeoff aborts, thus increasing engine operating costs. 
5 [016] In one embodiment, which can best be described with reference to a 
flow diagram 10 of Figures 1A-B, a core speed residual 11, an exhausted gas 
temperature residual 13, and a fuel flow residual 15 can be derived using an 
engine model, such as an empirical engine model or a first principle model. The 
first principle engine model may be based on a set of differential equations 

10 representing dynamics of the engine undergoing analysis and fault diagnosis. 
The empirical engine model may be expressed by a polynomial function of 
engine fan speed and ambient conditions such as ambient temperature, altitude 
and mach number. The engine model may be developed for a defined air bleed 
condition during takeoff, for example. An air bleed system removes 

15 compressed air from a compressor intermediate stage and bleeds the air into 
the atmosphere. During rapid transients, such as accelerations and 
decelerations, the bleed system is automatically opened and maintained open 
until the selected engine speed or power has been reached. 
[017] Alternatively, the empirical engine model may be expressed as a 

20 neural network model of the engine operating parameters. A neural network is 
an information processing paradigm similar to the method by which biological 
nervous systems, such as the brain, process information. Neural networks are 
characterized by powerful pattern matching and predictive capabilities in which 
input variables interact heavily to provide for modeling of nonlinear, complex 

25 dynamic systems. The neural network comprises a large number of highly 
interconnected processing elements (e.g., neurons) working in unison to solve 
specific problems. In biological systems the process of "learning" involves 
adjustments to the synaptic connections that exist between the neurons. In the 
present method, the neural network is not explicitly programmed but rather is 
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trained through exposure to real-time or to historical data. 
[018] A set of comparisons can be made, at constant engine fan speed 
and given ambient conditions, by taking measured vital engine performance 
parameters (i.e., engine operating data collected in the field by acquiring a core 
5 speed reading, an exhausted gas temperature reading, and a fuel flow reading) 
for an engine being analyzed, and comparing the readings with a set of 
corresponding standard engine characteristics or performance parameters for 
the engine model. The numerical differences between these two sets of 
performance parameters yield the core speed residual 11, the exhausted gas 

10 temperature residual 13, and the fuel flow residual 15. Recorded data 
comprising the engine performance parameters can be corrected to Sea Level 
International Standard Day conditions using altitude data so as to remove the 
effects of varying environmental conditions. Derivation of the residuals can be 
represented as, 

1 5 residual = (datum field - datum empirical ) 

[019] The mean of the core speed residual 11 and the standard deviation 
of the core speed residual 11 can be computed at step 17. The mean of the 
core speed residual 1 1 can be normalized to zero and the standard deviation of 
the core speed residual 11 can be normalized to unity, at step 19, to obtain a 

20 normalized core speed residual 21 by using normalization factors. 
Normalization factors which are equivalent to the mean and standard deviation 
obtained at step 17 can be obtained from the parameter distribution of a 
normally-operating engine, herein referred to as a baseline engine. 
[020] The normalization factors may be dynamically updated by obtaining 

25 engine performance parameters when the engine condition is determined to be 
normal. In accordance with the present method, the normalization factors are 
continually updated as additional engine performance parameters are acquired 
from the normally operating engine obtained in the field and during routine 
engine inspections. Engine performance parameters obtained from an engine 
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determined to be faulty are not included in an updated normalization factor. 
[021] An updated mean normalization factor may be obtained by: i) 
multiplying the current mean normalization factor by a first fraction to obtain a 
first product, ii) obtaining a current engine parameter from a normally-operating 
5 turbine engine, iii) multiplying the current engine parameter by a second fraction 
to obtain a second product, and iv) adding the first and second products to yield 
the updated mean deviation normalization factor, in accordance with the 
following expression, 

10 NFU~NF\y~x P 

where NF1 is a mean normalization factor, 

subscript p stands for a current time step, 

subscript p+1 stands for a step subsequent to the current time 
step p, 

15 n is the total number of samples accumulated at time step p, and 

x p are updated engine parameters. 
[022] An updated standard deviation normalization factor may be obtained 
from a current standard deviation normalization factor by: i) multiplying the 
square of the current standard deviation normalization factor by a third fraction 

20 to obtain a third product, ii) subtracting the updated mean normalization factor 
from the current engine parameter and multiplying the square of the result by a 
fourth fraction to obtain a fourth product, iii) subtracting the current mean 
normalization factor from the current engine parameter and multiplying the 
square of the result by the second fraction to obtain a fifth product, and iv) 

25 taking the square root of the sum of the third, fourth, and fifth products to yield 
the updated standard deviation, in accordance with the following expression, 
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NF2 P « = 



n-2 



NF2\ * ^ ix. ~ NF\J ♦ *iZM± 




n-l 



where NF2 is a standard deviation normalization factor. The mean and 
standard deviations may be used as slope and intercept coefficients, 
5 respectively, in a linear transformation, at step 19. 



store all the engine parameters to obtain an updated normalization factor. The 
only information required to update the normalization factors are the updated 
engine parameters and the normalization factors obtained in a previous time 
10 step. Accordingly, normalization factors can be updated dynamically resulting 
in a more accurate estimation of a baseline engine distribution without requiring 
large memory resources in which to store previously-obtained engine operating 
parameters. 

[024] As may further be appreciated by one skilled in the relevant art, the 
15 baseline engine distribution can differ from one engine to another because of 
variations in the respective engine component fabrication tolerances and in the 
manufacturing procedures for the engines. Therefore, an individual baseline 
engine distribution may be obtained for each turbine engine undergoing 
analysis. 

20 [025] The mean of the exhausted gas temperature residual 13 and the 
standard deviation of the exhausted gas temperature residual 13 can be 
computed at step 23. The mean of the exhausted gas temperature residual 13 
can be normalized to zero and the standard deviation of the exhausted gas 
temperature residual 13 can be normalized to unity, at step 25, to obtain a 

25 normalized exhaust gas temperature residual 27. The mean of the fuel flow 
residual 15 and the standard deviation of the fuel flow residual 15 can be 
computed at step 29. The mean of the fuel flow residual 15 can be normalized 
to zero and the standard deviation of the fuel flow residual 15 can be 



[023] 



As can be appreciated by one skilled in the art, it is not necessary to 
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normalized to unity, at step 31 , to obtain a normalized fuel flow residual 33. 
[026] The normalized core speed residual 21, the normalized exhausted 
gas temperature residual 27, and the normalized fuel flow residual 33 can be 
inputted as three corresponding input vectors to a clustering technique 
5 mapping, at step 35. As understood by one skilled in the art, a clustering 
technique mapping involves the task of grouping data points into homogeneous 
classes Or "clusters" so that data points in the same cluster are substantially 
similar and data points in different clusters are essentially dissimilar. 
[027] A primary concept in the clustering of data points is that a cluster 
10 comprises points separated by distances which are small in comparison to the 

distances between clusters. Alternatively, clustering can be thought of as a 

i 

form of data compression, where a large number of data points are mapped 
into a smaller number of representative clusters. The degree of similarity and 
of difference between individual observations are used to define the clusters, 

1 5 and to assign data points to clusters at varying degrees of aggregation. 

[028] The clustering technique mapping, in step 35, may comprise a self- 
organizing mapping operation. As understood by one skilled in the relevant art, 
a self-organizing map is a class of neural network which is 'trained 1 to structure 
itself using input signals. See, for example, Teuvo Kohonen, "The Self- 

20 Organizing Map," Proceedings of the IEEE, Volume 78, Number 9, September 
1990, pp. 1464-1480, incorporated herein by reference. 

[029] In the method of Figures 1A-B, the self-organizing map, which is a 
specific type of clustering algorithm forming neurons on a regular grid, may 
comprise a plurality of clusters, wherein each cluster may indicate or represent 
25 vector engine conditions, including both normal and faulty vector engine 
conditions. The self-organizing map 'learns* to recognize groups of similar input 
vectors in such a way that neurons physically near one another in a neuron 
layer respond to similar input vectors. 

[030] In an alternative embodiment, the clustering technique mapping may 
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comprise adaptive resonance theory, in which a central feature of an adaptive 
resonance theory system is a pattern-matching process, or vector classifier, 
that compares an external data input (i.e., an input vector) with an internal 
memory of an active code (i.e., a stored pattern). The input vector is classified 
5 into a category depending on the stored pattern the vector most closely 
resembles Adaptive resonance theory matching leads either to a resonant 
state or to a parallel memory search. If the search ends at an established 
code, a memory representation of the data may either remain the same or 
incorporate new data from matched portions of a current data input. If the 

10 search ends at a new code, the memory representation "learns" the current 
data input. In other words, once a stored pattern is found, the pattern is 
modified to resemble the input vector. If the input vector does not match any 
stored pattern within a certain predetermined tolerance, then a new category is 
created in memory by storing a new pattern similar to the original input vector. 

15 [031] In another alternate embodiment, the clustering technique mapping 
may comprise a least-squares partitioning method known as a K-means 
algorithm. The K-means algorithm functions by dividing a collection of objects 
into a plurality of K groups. The K-means algorithm iterates over two steps. 
First, the mean of each cluster is computed. The initial assignment of points to 

20 clusters can be done randomly. Second, the distance of each point from each 
cluster mean is computed, and each point is assigned to a nearest cluster. The 
iteration is repeated over the above two steps until the sum of "squared within 
group errors" (i.e., the distances of the points to the respective group means) is 
minimized and the function converges. The cluster groups obtained by this 

25 method are essentially as geometrically compact as possible around respective 
cluster means. 

[032] In yet another alternate embodiment, the clustering technique 
mapping may comprise a method employing a Gaussian mixture model. The 
Gaussian mixture method provides a probabilistic version of the traditional K- 
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means algorithm. However, whereas the traditional K-means algorithm 
represents clusters only in terms of the cluster centroids, the Gaussian mixture 
model operates on the basis that each cluster corresponds to an identifiable 
distribution of samples (i.e., engine parameters). Th Gaussian mixture method 
5 comprises a two-step iteration. First, via an "expectation step," an initial 
quantity of 'k' arbitrary centroids may be defined and a partial assignment of 
samples to clusters may be based on probabilities that the Gaussian mixture 
model assigns to the plurality of engine parameter samples. Second, via a 
"maximization step," the 'k' centroids are adjusted in conformance to the 

10 weighted means of the assigned samples. This two-step iteration is repeated 
until the sizes of adjustments required fall below a predefined threshold. 
[033] Accordingly, the clustering technique mapping, at step 35, can map 
an observation comprising a set of three input vectors (e.g., the normalized 
core speed residual 21, the normalized exhausted gas temperature residual 27, 

15 and the normalized fuel flow residual 33) into a nearest cluster to obtain a 
matching cluster number, at step 37. By interpreting the matching cluster 
number, the condition of the engine being analyzed can be determined and a 
failure (if present) can be diagnosed. The clustering technique mapping of 
residuals into engine condition space can be represented as, 

CTM 

20 normalized residual => engine condition space 

[034] The engine condition space may include various failure modes of 
interest including, for example, combustor liner burn-through failure, bleed band 
leakage failure, and exhaust gas temperature sensor failure. In one 
embodiment of the present invention, which used a self-organizing mapping 

25 method, the self-organizing map was trained for 1000 epochs, including data 
obtained from ten different reference turbine engines, and produced the cluster 
number relationship to vector engine condition shown in Table 1. A normal 
vector engine condition is indicated by a 'no trigger' condition. 
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Table 1 



Cluster Number 


Vector Engine Condition 


1 through 9 


No Trigger 


10 


EGT Sensor Failure Trigger (measuring 
lower than it should be) 


11 


EGT Sensor Failure Trigger (measuring 
higher than it should be) 


12 through 24 


No Trigger 


25 


Combustor Liner Failure Trigger 


26 


Bleed Band Failure Trigger 



[035] The combustor liner is a thermal coating disposed inside the engine 
combustor. The purpose of the combustor liner is to prevent burning fuel from 
5 overheating the walls of the combustion housing. A combustor-related failure 
may include cracking of the liner or a burn-through failure which could occur as 
a result of inferior baffle attachment welds or of a coked fuel nozzle resulting in 
poor fuel distribution and hot spots in a combustion chamber. 
[036] Bleed band related failures may generally occur more often than 

1 0 combustor liner failures. The most frequent failure mode of bleed band is bleed 
band leakage due to wear or improper sealing along the bleed band edge. 
Exhaust gas temperature sensors measure gas temperature at the exhaust of 
the turbine engine. In a typical sensor application, two thermocouple harness 
assemblies, each assembly having eight probes arranged on the circumference 

15 of a turbine nozzle containment ring, provide temperature signals to an engine 
electronic control unit. 

[037] At step 39, the matching cluster number may be used to determine 
whether the engine under analysis is normal or faulty. If the matching cluster 
number, found in Table 1, represents a faulty vector engine condition, this may 
20 be taken as an indication that the engine undergoing analysis may have a faulty 
condition and the diagnosis may conclude that the engine condition is 
determined to be faulty. If the matching cluster number represents a normal 
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vector engine condition the engine diagnosis may conclude that the engine 
condition is determined to be normal. The matching cluster number may also 
be used to compute a distance from the observation to the matching cluster 
number, at step 41. Using the information obtained in steps 39 and 41, a 
5 diagnosis confidence level, or belief factor, may be obtained, at step 43. That 
is, the closer the observation, or set of input vectors, is to the cluster having the 
matching cluster number, the higher the confidence level that the observation 
does belong to the group corresponding to the cluster having the matching 
cluster number. 

10 [038] The belief factor can be computed in either of two ways. If the 
matching cluster number indicates that the engine condition is faulty, at step 39, 
the belief factor may be obtained by subtracting from unity the ratio of i) the 
distance between the current observation and the closest cluster identified by 
the matching cluster number (i.e., a closest distance) and ii) the distance 

1 5 between the current observation and the next closest cluster (i.e., a next-closest 
distance). Thus, for a faulty engine condition determination, the belief factor 
can be given by, 



where 

20 x t is the observation, 

m c is the matching cluster closest to x i , and 

m nc is the next-closest cluster to x i 

[039] If the matching cluster number indicates that the engine condition is 
normal, at step 39, the belief factor may be obtained by subtracting from unity 
25 the largest ratio among the set of ratios obtained by dividing i) the closest 
distance (i.e., between the current observation and the matching cluster) by 
ii) the fault distances between the current observation and each of the clusters 
in the set of all the clusters representing a faulty condition. That is, for a normal 
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engine condition the following value is computed for a belief factor, 



BF. n =l-max r 



f*\ II ^ 



\x t -m 



f\ 



where m f is all clusters representing failures 



x f -m / | is the fault distance between the current observation and 

5 faulty condition cluster /. 

[040] As can be appreciated by one skilled in the relevant art, the 
disclosed method can be implemented by means of a computing module or 
device 50 such as shown in Fig. 2. The computing module 50 may comprise an 
application 51 in communication with a memory 53. The application 51 and the 

10 memory 53 may be provided on a computer drive (not shown) or may be 
provided via a computer-readable medium (not shown). The application 51 
may be executed by a processor 55. The memory 53 may include fault 
diagnosis data such as the engine operating data collected in the field, the 
normalization factors, the baseline engine parameter distribution, and the 

1 5 derived belief factor. Some or all of the fault diagnosis data may be transmitted 
to a remote storage device 57 via an interface 59, which can be a wired or 
wireless connection such as the Internet. Alternatively, the fault diagnosis data 
may be resident in the storage device 57 for access by the application 51 and 
processor 55 as required. 

20 [041] It should be understood, of course, that the foregoing relates to 
exemplary embodiments of the invention and that modifications may be made 
without departing from the spirit and scope of the invention as set forth in the 
following claims. 
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